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摘 要 : 半 高 宽 法 和 像 运动 法 是 估算 大 气 视 宁 度 的 常用 方法 。 这 两 种 方法 的 估算 精度 都 
会 受到 望远镜 跟踪 误差 ， 风 或 机 械 振动 等 因素 的 影响 。 提 高 视 宁 度 估算 精度 在 高 分 辨 成 像 、 
选 址 、 台 站 视 宁 度 监 测 等 领域 有 重要 意义 。 对 于 半 高 宽 法 ， 用 二 维 高 斯 函数 拟 合 长 曝光 点 源 
单 星星 像 并 选取 合适 方向 的 半 高 宽 估 算 视 宁 度 ， 与 谱 比 法 对 照发 现 可 以 有 效 改 善 跟踪 误差 ， 
风 或 机 械 振动 等 因素 的 影响 。 对 于 像 运动 法 ， 用 主 成 分 分 析 方法 找到 一 系列 斑点 图 重心 位 置 
在 合适 方向 上 投影 的 方差 ， 以 此 估算 视 宁 度 ， 与 谱 比 法 对 照发 现 可 以 有 效 改 善 跟 踪 误差 ， 风 
或 机 械 振 动 等 因素 的 影响 。 
关键 词 : MPA, 半 高 宽 ; 像 运动 法 ; 跟踪 误差 


中 图 分 类 号 ; P11 文献 标识 码 :A 
1 引 E 


大 气 满 流 是 降低 地 面 大 型 光学 望远镜 成 像 质 量 的 主要 因素 之 一 , 它 会 使 图 像 在 像 面 发 生 
随机 的 拌 动 ， 畸 变 和 破碎 。 为 描述 大 气 满 流 对 成 像 过 程 的 影响 ， 基 于 Komolgorov 关于 大 气 
注 流 的 统计 模型 ，D.L.Friedm1 引 入 了 视 宁 度 参数 mm， 又 被 称 为 Fried 参数 ， 或 者 大 气相 干 长 
度 。 根 据 定义 ，ro 表 征 着 地 面 望远镜 实现 衍射 极限 成 像 的 极限 口径 ">， 也 代表 着 大 气 淇 流 的 
剧烈 程度 。ro 值 在 选 址 和 高 分 辨 成 像 等 领域 有 极 高 的 重要 性 ， 比 如 斑点 干涉 术 “ 的 理论 传递 
函数 模型 与 yp 值 有 关 "， 自 适应 光学 中 Shack-Hartmann 波 前 传感器 微 透镜 阵列 的 数量 由 ro 
所 决定 ”。 
HEMT ro 有 许多 种 方法 。 有 些 需 要 使 用 专门 的 仪器 ， 比 如 ， 利 用 北极 星星 轨 的 横 
向 拌 动 方差 估算 视 宁 度 ， 利 用 探 空气 球 或 者 声 雷 达 测 量 不 同 高 度 大 气 的 温度 结构 常数 从 
而 推导 出 大 气 视 宁 度 ” "， 利 用 刀口 前 后 焦 面 光 强 变 化 程度 从 而 推导 出 视 宁 度 ，”， 利 用 干涉 
仪 中 被 大 气 干 扰 的 星 像 的 干涉 条 纹 来 推导 出 视 宁 度 ”, 利用 望远镜 上 两 个 间距 超过 ro FFL 
径 的 差分 到 达 角 方差 来 推算 视 宁 度 汪 ， 利 用 星 像 闪烁 的 程度 来 推算 视 宁 度 “ “等 等 。 

有 些 则 可 以 利用 望远镜 观测 过 程 中 的 次 点 图 (曝光 时 间 在 大 气相 干 时 间 量 级 的 短 
上 曝光 图 像 ) 估 算 ， 比 如 ， 将 一 组 点 源 星 像 的 实际 平均 长 曝光 或 者 平均 短 曝 光 传递 函数 "与 对 
应 的 理论 传递 函数 相 比较 从 而 推导 出 视 宁 度 ” , 通过 将 一 组 图 像 的 实际 谱 比值 与 理论 值 相 比 
较 来 得 到 视 宁 度 ” (在 扩展 源 图 像 的 高 分 辩 重 建 领域 应 用 广泛 , 本文 称 为 谱 比 法 ) , 利用 大 气 
注 流 对 图 像 对 比 度 造成 的 变化 从 而 反 推出 视 宁 度 “， 直 接 利 点 源 目标 长 曝光 星 像 的 半 高 宽 
(FWHM) 估算 视 宁 度 “”"”，, 利用 星 像 在 像 面 抖动 的 方差 来 估算 视 宁 度 ” (本 文 称 为 像 运 动 法 
(ImageMotion) ) 等 等 。 

在 诸多 方法 中 ，FWHM 法 和 像 运 动 法 可 以 较为 方便 快捷 的 估算 出 视 宁 度 参 数 ， 是 应 用 
十 分 广泛 的 方法 。 当 望远镜 受到 跟踪 误差 ， 风 或 者 机 械 振动 等 因素 (以 下 统称 为 综合 因素 ) 
的 影响 时 ， 这 两 种 方法 往往 会 受到 影响 ， 从 而 降低 视 宁 度 的 估算 精度 。 
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本 文 将 介绍 FWHM 法 估算 视 守 度 的 基本 原理 ,分 析 其 应 用 条 件 ， 及 综合 因素 将 产生 何 
种 影响 ,并 给 出 应 对 方法 。 介 绍 像 运 动 法 的 基本 原理 ,跟踪 误差 等 因素 的 影响 及 给 出 应 对 方 
法 。 随 后 将 应 对 方法 应 用 到 实际 数据 中 去 ， 并 分 析 取 得 的 结果 ， 将 结果 与 用 谱 比 法 对 照 ， 以 
说 明 应 对 方法 的 有 效 性 。 


2 方 法 


在 等 晕 区 内 ， 大 气 - 望 远 镜 综合 系统 可 以 看 做 一 个 线性 空 不 变 系统 ， 短 曝光 成 像 过 程 在 
空域 和 频 域 可 分 别 用 (1) 式 和 (2) 式 表示 表示 : 
i(x) = o(x) * s(x) 
1(q) = O(q)S(q) 
式 中 ，;o) 表 示 短 曝光 图 像 ，1fg) 是 其 对 应 的 频 域 表 示 ，of(xz) 表 示 目 标 ，Of(g) 表 示 其 傅 里 时 变 


换 ，s(x) 表 示 点 扩散 函数 ， 对 于 点 源 目标 来 说 i(x)=s(x)，S(q) 表 示 短 曝光 大 气 - 望 远 镜 综合 系 


统 的 传递 函数 ，g 表示 归 一 化 空间 频率 ，* 表 示 卷 积 。 
短 曝 光 大气- 望远镜 综合 系统 的 传递 函数 又 可 以 分 解 为 望远镜 的 传递 函数 与 短 曝 光大 气 
的 传递 函数 之 积 : 


S(q) =T(q)B(q) 
2.1 i 9E (FWHM) 法 
半 高 宽 法 的 本 质 是 利用 点 源 单 星 长 曝光 图 像 的 半 高 宽 与 视 宁 度 参数 之 间 的 简单 关系 : 
FWHM = 0.934 


h 


(4) 
其 中 4 表示 波长 。 长 曝光 图 像 可 以 通过 多 张 短 上 曝光 图 像 的 统计 平均 来 获得 ”。 因 为 长 曝光 星 
像 的 强度 分 布 近似 于 高 斯 函数 ”， 常 常 使 用 高 斯 函数 拟 合 以 计算 半 高 宽 ,高 斯 函数 的 半 
高 宽 与 标准 差 存 在 一 简单 关系 : 


FWHM = 2V21n20 


分 析 长 曝光 大 气 - 望 远 镜 综合 系统 的 分 辩 率 与 大 气 的 分 辨 率 的 比值 随 D/ ro 的 变化 情况 ， 
分 辨 率 定 义 为 在 频率 域 对 MTF 积分 ”””。 发 现 当 望远镜 口径 D 远大 于 视 宁 度 ww 时 ， 系 统 
分 辨 率 越 来 越 接 近 大 气 的 分 辩 率 ， 因 此 ， 当 D > no 时 , 可 以 将 长 曝光 综合 系统 的 传递 函数 与 
长 曝光 大 气 的 传递 函数 划 等 号 。 


(SD); * (BD), =exp(-3.44{ L] 


()RAR REE, a-ryD , LE 表示 长 曝光 。 


为 了 获取 空域 长 曝光 点 扩散 函数 (so0)7z 的 半 高 宽 ， 需 将 (6) 式 做 反 傅 里 叶 变 换 ， 再 取 其 
半 高 宽 ， 但 是 ，5/3 次 方 使 我 们 难以 得 到 解析 解 ， 一 种 近似 方法 是 用 2 次 来 替代 5/3 次 以 用 
高 斯 函数 来 近似 传递 函数 ， 如 此 即 可 得 到 半 高 宽 的 解析 式 (4) 。 

根据 上 述 分 析 ， 公 式 (4) 的 推导 过 程 实际 上 做 了 两 处 近似 ，1. 用 大 气 的 传递 函数 代替 系 
统 的 传递 函数 ，2. 用 2 次 方 替 代 5/3 次 方 。 为 了 说 明 两 处 近似 在 何 种 程度 上 影响 着 (4) 式 的 
准确 度 ， 使 用 数值 方法 计算 出 理论 上 系统 长 曝光 点 扩散 函数 的 半 高 宽 然 后 将 其 与 (4) 式 的 半 
高 宽 做 比较 。 其 结果 如 图 1 所 示 。 
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图 1 (4) 式 半 高 宽 与 综合 系统 的 半 高 宽 的 比值 随 D/ro 的 变化 情况 
Fig.1 Ratio of FWHM calculated by equation (4) and theoretical FWHM of the synthetic system as a function of D/ ro 
1 说 明 ， 当 D > ro, (4) 式 才能 达到 较 高 的 准确 度 ， 比 如 当 D=Im, ro=10cm 时 ， 
半 高 宽 比 值 可 以 达到 0. 93。 说 明 对 于 口径 Im 级 别 的 望远镜 ，(4) 式 在 理论 上 已 经 能 够 达到 


较 高 的 准确 度 ， 因 此 本 文 使 用 (4) 式 估算 roo 


上 述 分 析 是 基于 假设 望远镜 是 理想 望远镜 的 情况 ， 即 无 像 差 , 无 跟踪 误差 , JER 
响 ， 支 撑 系 统 完全 刚性 等 。 实际 上 ， 这 些 因素 痢 会 使 得 实际 拍摄 的 长 曝光 点 源 星 像 偏 大 ， 从 
而 使 得 公式 (4) 估算 的 wo 值 偏 小 。 

当 受 到 综合 因素 影响 时 , 长 曝光 图 像 往往 会 在 某 个 方向 上 延展 开 来 ,强度 分 布 表现 出 类 
椭圆 形 ， 选择 合适 的 方向 来 计算 半 高 宽 变 得 尤为 重要 。 应 对 方法 是 使 用 二 维 高 斯 函数 拟 合 受 
到 综合 因素 影响 的 长 曝光 点 源 单 星 图 像 ,将 方差 最 大 的 方向 视 为 受到 影响 较 大 的 方向 。 用 于 
拟 合 的 二 维 高 斯 函数 如 下 所 示 : 


f (al, x2) = Aexp( 


1 2 1 £4 2 F2 
Ce 27 quy aper (aaa) 
20d-p 2) O, (en O, 


o; poo, 
90,0, 0; 
绝 大 多 数 情况 下 ， 相 关系 数 o NAS, WIENER FA EME, ThE) 25 BE 
对 角 和 矩阵 以 找到 受到 影响 较 大 的 方向 ,再 取 与 之 垂直 方向 上 的 方差 , 即 对 角 和 矩阵 中 较 小 的 方 


对 应 的 协 方差 矩阵 为 : 


差 ， 以 此 来 计算 半 高 宽 ， 进 而 使 用 公式 (4) 估算 ro WE 2 所 示 。 
1754 175 
150 150 
125 125 

Om 
F 1 ES F 1 
n 75 ý 75 
504 504 On 

25 25 

9 十 一 0 十 一 一 —— 
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图 2 左 图 是 一 张 受 到 综合 因素 影响 的 点 源 单 星 的 长 曝光 图 像 ， 右 图 是 拟 合 长 曝光 图 像 的 二 维 高 斯 函数 图 像 


Fig.2 A long exposure image affected by tracking error, wind or mechanical vibration(Left), A two dimensional gaussian function 
that fits the left long exposure image(Right) 
2.2 像 运动 法 (ImageMotion) 
大 气 满 流 对 成 像 过 程 的 一 种 影响 就 是 改变 瞳 面 的 波 前 到 达 角 , 造成 像 面 图 像 的 随机 抖动 ， 
因此 测量 图 像 拌 动 的 剧烈 程度 就 可 以 间接 地 推导 出 视 宁 度 参数 roo D.L.Fried 揭示 了 ro 与 
波 前 到 达 角 方差 之 间 的 关系 : 


? 9.458(2)3(2)3 
on - 3583 CD 
FEA D RANIG LAAN HUS. o, 表示 波 前 到 达 角 标准 差 。 
对 于 点 源 单 星 图 像 , 重心 的 位 置 变化 就 对 应 着 瞳 面 的 波 前 到 达 角 起 伏 ”。 因 此 可 以 根据 
计算 像 面 图 像 重 心 的 抖动 方差 推算 出 视 宁 度 参数 ry。 其 中 ， 重 心 的 计算 方法 如 下 ™: 


二 


人 (OM ») 


其 中 xy RARE OLE, (xy) AN Go y) A b ABB ER Ho X 03 Y EOD be 

像 运动 法 会 受到 综合 因素 的 影响 。 假 设 仅 存在 大 气 清流 影响 时 图 像 重 心 位 置 为 Xa Ek 
合 因素 造成 的 偏 移 量 为 Xy， 有 是 X4 与 XYy 相 互 独立 ， 则 实际 重心 位 置 为 X=X4+Xy， 易 得 出 实 
际 重 心 位 置 的 方差 DX=DX,+DXy, W) DX 宇 DX4。 因 此 ， 一 般 来 说 ， 综 合 因素 会 造成 图 像 拌 
动 方差 增 大 ， 从 而 使 得 估算 的 ro Emh o 

每 张 图 的 曝光 时 间 长 度 也 会 影响 到 像 运 动 法 ,理想 情况 下 ,曝光 时 长 小 于 或 等 于 大 
气相 干 时 间 ， 不 同 斑 点 图 对 应 着 不 同 状 态 的 “冻结 ”的 大 气 ， 图 像 重 心 位 置 也 不 一 样 。 当 每 
张 图 的 曝光 时 长 增加 时 ,重心 位 置 会 被 平滑 掉 ， 从 而 降低 重心 位 置 拌 动 的 方差 ， 造成 估算 的 
ro (MWK © 

需要 指出 的 是 ， 望 远 镜 口径 大 小 也 会 对 像 运 动 法 造成 影响 。H. M. Martinpg 认 为 像 运 动 
实际 上 是 由 大 于 口径 的 汕 流 元 胞 造成 的 , 小 尺度 的 汕 流 尽管 会 造成 图 像 的 斑点 结构 ,但 是 口 
径 内 众多 的 小 尺度 满 流 会 平滑 掉 重 心 位 置 , 造成 重心 位 置 偏 移 不 明显 。G. Ricort! 1 使 用 40cm 
口径 的 望远镜 估算 视 守 度 ， 结 果 显 示 像 运动 法 估算 的 视 宁 度 是 对 比 度 法 的 1.4 倍 。P R. 
GOODE” (FH 65cm 口径 的 望远镜 估算 视 宁 度 ， 结 果 显 示 像 运动 法 的 结果 是 谱 比 法 的 1.5 
fit. Bin Ma 使 用 50cm 口径 的 望远镜 对 比 了 星 轨 法 (本 质 上 是 像 运动 法 ) 与 DIMM， 发 现 星 
轨 法 的 结果 偏 大 。 这 些 研究 表明 ， 大 口径 的 望远镜 使 用 像 运 动 法 估算 视 宁 度 , 会 导致 结果 偏 
大 。 

公式 (9) 是 根据 kolmogorov 谱 推 导 的 , 并 没有 考虑 大 气 灌流 外 尺度 的 影响 , 实际 上 , Von 
Karman 谱 的 数值 模拟 表明 ， 当 考虑 有 限 外 尺度 时 ， 重 心 的 抖动 方差 会 偏 小 ， 导 致 估算 的 视 
宁 度 偏 大 。 

当 图 像 信 品 比较 低 时 ，CCD 读 出 噪声 会 影响 图 像 重 心 位 置 的 计算 ， 造 成 像 运 动 法 视 宁 
度 的 不 准确 估算 ， 通 常 的 应 对 方法 有 : 1. 门限 值 法 ， 对 图 像 设置 门限 值 ( 一 般 使 用 图 像 最 大 
值 乘 以 某 个 小 于 1 的 系数 来 表示 ) ， 仅 仅 使 用 门限 值 以 上 的 像素 值 计算 重 心 ，2. 加 窗 法 ， 设 
置 一 个 窗口 范围 , 仅仅 使 用 窗口 内 部 的 像素 值 计 算 重 心 。 如 果 对 一 组 图 像 直接 使 用 门限 值 法 
估算 视 宁 度 ， 可 以 发 现 门限 值 的 大 小 会 严重 影响 到 估算 结果 。 如 图 3 所 示 : 


im] 


(9) 


(10) 
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图 3 门限 值 大 小 对 视 宁 度 估算 的 影响 
Fig .3 Influence of threshold on seeing estimation 

J. Verin 等 人 分 析 了 CCD 噪声 对 重心 位 置 影响 ， 发 现 重心 位 置 抖动 方差 与 窗口 大 小 
呈 四 次 方 关系 ,是 严重 影响 因素 ， 设 置 门限 值 在 某 种 意义 上 就 是 限制 了 窗口 大 小 ,加 窗 法 则 
是 直接 就 将 窗口 限制 为 某 个 较 小 的 范围 。 

本 文 将 加 窗 法 与 门限 值 法 相 结 合 ， 首 先 以 图 像 最 大 值 为 中 心 ， 以 Wp AAI, Wo 
定义 为 使 用 高 斯 函数 近似 点 扩散 函数 之 后 ， 强 度 值 下 降 为 最 大 值 的 Le 时 的 直径 ， 类 似 于 衍 
射 受 限 点 扩散 函数 的 rms. 宽度 的 定义 ”。 理论 上 , Wp*12FWHMj- 随后 , 考虑 到 Sandrine 
Thomas"?.5 Chao Li 等 人 模拟 分 析 发 现 以 3 倍 的 读 出 噪声 作为 门限 值 可 以 有 效 地 降低 读 
出 噪声 的 影响 , 本 文 也 以 3 倍 的 读 出 噪声 作为 门限 值 计算 重心 位 置 , 读 出 噪声 可 以 从 一 系列 
上 暗 场 图 像 中 分 析 获 得 。 

当 受 到 综合 因素 的 影响 时 ,一 系列 点 源 单 星 斑 点 图 的 重心 位 置 往往 会 在 某 个 方向 延展 开 
来 ， 使 得 这 个 方向 上 的 抖动 方差 要 大 于 其 他 方向 。 应 对 策略 是 使 用 主 成 分 分 析 方 法 (PCA) 
找到 一 组 重心 坐标 数据 中 方差 最 大 的 方向 , 即 第 一 主 成 分 方向 , 将 此 方向 看 做 受到 影响 较 大 
的 方向 。 将 坐标 数据 降 维 在 第 二 主 成 分 方向 上 (与 第 一 主 成 分 正 交 的 方向 ) ， 并 计算 其 方差 ， 
然后 利用 该 方差 来 估算 wo 值 。 在 图 4 的 情况 下 ， 就 是 取 y 方 和 同上 的 方差 来 估算 roo 

主 成 分 分 析 方 法 是 一 种 应 用 广泛 的 数据 降 维 手段 , 在 天 文 数据 处 理 领 域 也 有 应 用 ”, 一 
般 是 将 数据 降 维 到 前 n 阶 方差 最 大 的 方向 上 以 保留 原始 数据 的 主要 信息 特征 , 本文 的 用 法 则 
是 将 二 维 的 重心 位 置 坐标 数据 降 维 到 方差 小 的 方向 上 。 具 体 的 算法 实现 如 下 : 

CD 将 所 有 重心 坐标 组 成 矩阵 


X=(' xj (11) 
» » Yn 
(2) 将 和 每 一 行进 行 零 均 值 化 
X,,-X-X (12) 
(3) 计算 出 协 方差 矩阵 
1 
=— Xx xl 1 
C N new new ( 3 ) 


(4) 将 协 方差 矩阵 C 对 角 化 为 C" 
(5) C 中 对 角 线 的 元 素 即 为 第 一 主 成 分 和 第 二 主 成 分 方向 的 方差 ， 值 小 的 方差 即 为 所 需 
要 的 方差 


202009.00006v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


-30 -20 -10 0 10 20 30 


图 4 一 系列 斑点 图 重心 坐标 位 置 及 第 一 主 成 分 方向 和 第 二 主 成 分 方向 


Fig.4 Coordinates of center of gravity of a series of speckle images and directions of the first and second principle component 

尽管 选取 了 方差 小 的 方向 来 估算 视 宁 度 , 但 是 并 非 认为 此 方向 完全 不 受到 任何 影响 , 有 
必要 限制 一 组 斑点 图 的 拍摄 时 间 ， 因 为 长 时 间 风 向 等 因素 可 能 发 生 剧烈 变化 ,导致 各 个 方向 
上 都 会 受到 影响 。 可 以 用 以 下 方法 判断 所 选择 的 方向 是 否 也 受到 了 较 大 的 干扰 : 将 重心 坐标 
在 第 二 主 成 分 方向 的 投影 以 直方 图 的 形式 表示 ， 并 使 用 高 斯 函数 拟 合 ， 并 考察 其 决定 系数 ， 
因为 当 不 存在 任何 影响 时 , 重心 位 置 的 变化 应 当 符 合 高 斯 分 布 , 受到 影响 较 小 的 方向 也 应 该 
符合 高 斯 分 布 。 实 际 应 用 中 可 以 限制 决定 系数 ， 比 如 当 决 定 系 数 大 于 0.8 时 才 使 用 本 文 的 方 
法 ， 和 否则 应 当 放弃 使 用 该 组 图 像 估算 。 图 5 为 一 组 斑点 图 在 第 二 主 成 分 上 投影 的 直方 图 。 


projected center of gravity histogram 


0 5 
R? =0.925119 


5 在 第 二 主 成 分 方向 上 投影 的 重心 坐标 以 及 拟 合 的 高 斯 函数 
Fig .5 Coordinates of center of gravity projected on the second PCA direction and fiited gaussian function 
2.3 谱 比 法 
谱 比 法 由 0. Von Der Lühe 提出 ”"， 是 通过 计算 一 组 图 像 平均 频谱 的 平方 与 平均 能 谱 的 
比值 ， 再 将 此 值 与 理论 值 作 比 较 ， 对 于 给 定 的 望远镜 来 说 ， 理 论 值 仅仅 是 视 宁 度 参数 ro 的 
函数 。 其 表达 式 如 下 : 


ge. M) ZOD _(S@y 
(^0) Pola) (S) 
HP, (S(g) BERT ERST 3 TONG 38 en Hc lo T ECT f OCHO ERG, EDT TOT 
传递 函数 时 , 意味 着 将 所 有 图 像 按照 重心 对 齐 再 做 计算 , 因此 可 以 不 受 跟踪 误差 等 因素 的 影 
响 ， 表 示 平 均 短 曝光 传递 函数 时 ， 其 理论 表达 式 如 下 : 


(14) 


Hy 


(S(q))., -T'q)exp(-688 ) d -q*) 


平均 能 谱 的 理论 表达 式 如 下 : 


(s*@)=—Jou.ms.ndy 


pus 
3 


QCy.q) = e068 [a aa +(y-a)3 - (y *q|- à» * Qv - a]: op 


Sog) = [wc 


N Kae 化 常数 ， 与 4 都 是 模 的 范 


人 比值 是 ro "n " 


线 比较 就 可 以 佑 算出 视 宁 度 参数 roo 


3 数据 处 理 与 分 析 


224 = 
Lor. pwe = DWG 


围 为 0 到 1 的 矢量 


镜 ，(S2(@) 只 与 WAR: ( 
pe ro 值 就 可 以 计算 理论 .| 


上 对 应 的 谱 比 曲线 ， 将 其 


Zo S(y,q) HUAFA 1 的 
S(r BRS 3 MAK; 


圆 的 交 


因此 R 的 
与 实 


际 谱 比 


3.1 数据 

为 验证 本 文 应 对 综合 因素 影响 的 方法 的 有 效 性 ， 我 们 采用 二 组 实测 数据 进行 试验 。 

第 一 组 数据 是 2014 年 9 月 12 日 晚 UT 18:48:28 到 UT 19:12:03 时 间 段 内 ， 抚 仙 湖 1 米 
新 真空 太阳 望远镜 (NVST) 所 拍摄 的 25 组 数据 ， 每 组 数据 含 200 张 斑 点 图 ， 详 细 信息 如 表 1 
所 示 : 

表 1 实验 参数 
Table 1 Experiment Parameters 
Parameter Value 
Aperture Diameter 0.98m 
Central Wavelength 705.8nm 
Pixel Size 6um/0.04 arcsec 
Target HIP9487 
Region of Interest Size 192x192 pixels 
Exposure Time 30ms 
Duration of Observation 23 min 34 sec 
Number of Groups 25 
Frames per Group 200 
Duration per Group 20 sec 
Zenith Angle 23° 4' 59” --25° 33 SS" 

第 二 组 数据 是 2009 年 2 月 9 日 晚 UT 12:49 到 UT 20:46 内 ， 使 用 丽江 高 美 古 观测 站 的 

2.4 米 望远镜 针对 10 个 目标 拍摄 的 95 组 图 像 ， 每 组 500 帧 。 详 细 信息 如 表 2 所 示 。 
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表 2 实验 参数 


Table 2 Experiment Parameters 


Parameter Value 
Aperture Diameter 2.4m 
Central Wavelength 550nm 
Pixel Size 0.021 arcsec 
Target HIP30960, HIP25245, HIP32993, HIP43635, HIP48670, 


HIP53360, HIP64543, HR4300, HR4335, HR4518 


Region of Interest Size 256x256 pixels 
Exposure Time 8ms 
Duration of Observation "]hours57minutes 
Number of Groups 95 
Frames per Group 500 
Duration per Group «10 sec 
Zenith Angle 8° 16' 12" —-21* 4' 3" 


对 所 有 图 像 按照 文献 ”中 的 方法 进行 预 处 理 。 
对 于 NVST 的 每 一 组 数据 ， 以 任意 一 张 图 (本文 选取 第 一 张 图 ) 的 最 大 值 为 中 心 ， 将 图 像 
裁剪 成 192x192 像素 的 小 子 图 ， 目 的 是 为 了 避免 非 等 晕 效 应 的 影响 和 加 快 处 理 速 度 。 
对 于 高 美 古 2. 4m 望远镜 的 数据 ， 因 为 图 像 数 据 已 经 是 以 5. 376x5. 376 角 秒 的 小 子 图 读 
出 的 ， 因 此 不 用 再 做 裁 切 处 理 了 。 

图 6 是 经 过 预 处 理 所 得 到 的 斑点 图 。 
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图 6 左边 为 校准 完成 的 NVST 拍摄 的 斑点 图 ， 右 边 为 校准 完成 的 2. Am 望远镜 拍摄 的 斑点 图 
Fig.6 a frame of calibrated speckle image taken by NVST(Left), a frame of calibrated speckle image taken by 2.4m tele- 
scope(Right) 

3.2 处 理 与 分 析 

对 所 获 数据 采用 以 上 三 种 方法 进行 处 理 。 详 细 处 理 流程 如 图 7 所 示 。 对 半 高 宽 法 ， 预 处 
理 完成 后 ,直接 将 所 有 斑点 图 全 加 ， 等 效 为 单 张 长 曝光 图 像 ， 再 对 该 图 像 使 用 二 维 高 斯 函数 
拟 合 , 并 使 用 上 述 方法 找到 合适 方向 的 半 高 宽 从 而 估算 视 宁 度 。 对 像 运动 法 , 预 处 理 完 成 后 ， 
计算 所 有 图 像 重心 坐标 ,并 使 用 第 二 节 所 述 方法 将 重心 坐标 投影 在 第 二 主 成 分 上 ， 以 直方 图 
表示 该 方向 上 重心 坐标 ， 计 算 决定 系数 ， 如 果 决 定 系 数 R 大 于 0.8， 就 使 用 该 方向 上 方差 
估算 视 宁 度 ， 和 否则 就 不 使 用 这 组 图 像 ， 这 个 条 件 也 应 用 于 半 高 宽 法 。 对 于 谱 比 法 ， 预 处 理 完 
成 后 , 将 所 有 点 源 单 星 图 像 按照 重心 对 齐 ， 随 后 以 对 齐 后 的 图 像 计 算 平 均 短 曝光 频谱 与 平均 
能 谱 ， 得 到 实际 谱 比 值 ， 与 理论 谱 比值 比较 ， 找 到 最 为 接近 的 视 宁 度 参 数值 。 
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为 了 方便 比较 , 将 所 估算 的 ro 全 部 归 算 到 天 顶 距 为 0”， 且 波长 为 500nm 时 的 值 。 归 算 


所 依据 的 公式 如 下 : 
A cosy 3 
[c aman 


中 ，4 为 波长 ，Y 为 天 项 距 ，C# 为 折射 率 结构 常数 。 
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7 三 种 方法 的 流程 图 


Fig .7 Flowchart of three methods 
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2. 


作为 对 比 , 同时 使 用 以 上 应 对 综合 因素 的 修正 方法 和 选取 任 一 方向 的 方差 来 估算 视 宁 度 


(本 文选 取 纵 向 的 方差 ) 。 对 于 NVST 的 数据 ， 使 用 纵向 方差 估算 的 视 宁 度 如 图 8 所 示 。 使 有 


H 


修正 方法 估算 的 视 宁 度 如 图 9 所 示 。 对 于 高 美 古 2. 4 望远镜 的 数据 ， 使 用 纵向 方差 估算 的 


视 宁 度 如 图 10 所 示 。 使 用 修正 方法 估算 的 视 宁 度 如 图 11 所 示 。 
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图 8 NVST 数据 ; 点 划 线 是 在 FWHM 法 中 使 用 纵向 半 高 宽 估算 的 视 宁 度 ， 虚 线 是 在 像 运 动 法 


度 ， 实 线 是 使 用 谱 比 法 估算 的 视 宁 度 


P fir 


纵向 方差 估算 的 视 宁 


Fig.8 NVST data: Dash-dotted line shows seeing parameter estimated by using longitudinal FWHM in the FWHM method, dotted 


line shows seeing parameter estimated by employing longitudinal variance in the Image Motion method, solid line shows seeing parame- 


ter estimated by spectral ratio technic. 
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图 9 NVST 数据 ; 点 划 线 是 使 用 修正 的 FWHM 法 估算 的 视 宁 度 ， 虚 线 是 使 用 修正 的 像 运 动 法 估算 的 视 宁 度 ， 实 线 是 使 


谱 比 法 估算 的 视 宁 度 


io 


Fig.9 NVST data: Dash-dotted line shows seeing parameter estimated using corrected FWHM method, dotted line shows seeing 


parameter estimated by corrected Image Motion method, solid line shows seeing parameter estimated by spectral ratio technic 
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图 10 2.4m 望远镜 数据 ; 点 划 线 是 在 FWHM 法 中 使 用 纵向 半 高 宽 估算 的 视 宁 度 ， 虚 线 是 在 像 运 动 法 中 使 用 纵向 方差 估算 
的 视 宁 度 ， 实 线 是 使 用 谱 比 法 估算 的 视 宁 度 


Fig.10 2.4m telescope data: Dash-dotted line shows seeing parameter estimated by using longitudinal FWHM in the FWHM 


method, dotted line shows seeing parameter estimated by employing longitudinal variance in the Image Motion method, solid line shows 


seeing parameter estimated by spectral ratio technic. 
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Group 


图 11 2.4m 望远镜 数据 ， 点 划 线 是 使 用 修正 的 FWHM 法 估算 的 视 宁 度 ， 虚 线 是 使 用 修正 的 像 运动 法 估算 的 视 宁 度 ， 实 线 
是 使 用 谱 比 法 估算 的 视 宁 度 


Fig.11 2.4m telescope data: Dash-dotted line shows seeing parameter estimated using corrected FWHM method, dotted line shows 


Seeing parameter estimated by corrected Image Motion method, solid line shows seeing parameter estimated by Spectral ratio technic 
图 9 显示 ,对 于 NVST 的 数据 ,在 观测 时 间 段 内 ,使 用 FWHM 法 得 到 的 加 均值 为 9. 96cm, 
标准 差 为 0. 73cm， 像 运动 法 得 到 的 ro 均值 为 17. 44cm， 标 准 差 为 2. 12cm， 谱 比 法 得 到 的 
ro 均值 为 11. 04cm， 标 准 差 为 0. 62cm。 图 11 显示 ， 对 于 高 美 古 2. 4 望远镜 的 数据 ， 使 用 
FWHM 法 得 到 的 ro 均值 为 7. 13cm， 标 准 差 为 0. 82cm， 像 运动 法 得 到 的 r0 均值 为 13. 48cm, 
标准 差 为 2cm， 谱 比 法 得 到 的 r0 均值 为 7. 88cm， 标 准 差 为 0. 97cm. 
从 以 上 四 幅 图 可 以 得 到 以 下 信息 : 
第 一 , FWHM 法 得 到 的 视 宁 度 比 谱 比 法 偏 小 , 可 能 的 原因 有 1. (4) 式 的 推导 过 程 本 身 就 说 


Chinay nj 人 人 天 甘 日 工 || 
CninaX iv& FEKT! 


明 ， 估 算 的 r0 值 会 稍微 偏 小 ，2. 望远镜 像 差 ， 残 余 的 跟踪 误差 ， 风 或 机 械 振 动因 素 仍 然 使 
得 半 高 宽 偏 大 ， 从 而 使 估算 出 的 ro 值 更 偏 小 。 

第 二 ， 像 运动 法 估算 的 视 宁 度 比 谱 比 法 偏 大 ， 可 能 的 原因 有 1. 曝光 时 间 过 长 ， 导 致 实 
际 测量 的 像 运动 方差 偏 小 ,造成 估算 的 视 宁 度 偏 大 ，2. 望远镜 口径 过 大 ， 像 运动 更 多 的 反应 
的 是 大 于 口径 的 满 流 元 胞 的 起 伏 , 3. 没有 考虑 有 限 外 尺度 的 影响 4. CCD 噪声 和 光子 噪声 对 重 
心 拌 动 方差 计算 的 影响 仍然 存在 。 
第 三 ， 纵 向 方差 估算 的 r0 比 修正 方法 估算 的 r0 偏 小 , 这 是 因为 纵向 方差 更 多 的 包括 了 
跟踪 误差 ， 风 或 机 械 振动 等 因素 的 影响 而 偏 大 ， 使 得 估算 的 r0 值 偏 小 
第 四 ,对 于 nvst 的 数据 ,修正 的 FWHM 法 估算 的 视 宁 度 与 谱 比 法 之 间 的 相关 系数 为 0. 79， 
而 纵向 半 高 宽 估算 的 视 宁 度 与 谱 比 法 之 间 的 相关 系数 为 0.39。 修 正 的 像 运 动 法 与 谱 比 法 之 
间 的 相关 系数 为 0.52， 而 纵向 方差 估算 的 视 宁 度 与 谱 比 法 之 间 的 相关 系数 为 -0.01。 对 于 高 
JET 2. 4 望远镜 的 数据 ， 修 正 的 FWHM 法 估算 的 视 宁 度 与 谱 比 法 之 间 的 相关 系数 为 0. 91， 
而 纵向 半 高 宽 估算 的 视 宁 度 与 谱 比 法 之 间 的 相关 系数 为 0. 58, o 修正 像 运动 法 与 谱 比 法 之 间 
的 相关 系数 为 0.65， 而 纵向 方差 估算 的 视 宁 度 之 间 的 相关 系数 为 0. 24。 这 说 明 ， 在 一 定时 
HA, 当 一 系列 斑点 图 受到 了 综合 因素 的 影响 , 而 本 文 所 选取 的 方向 并 未 受到 严重 影响 时 (最 
> 低 要 求 决 定 系 数 大 于 0. 8) ， 使 用 本 文 的 方法 可 以 有 效 提高 半 高 宽 法 ， 像 运动 法 与 谱 比 法 之 
闻 的 相关 性 ,对 比 任意 选择 一 个 方向 估算 视 宁 度 , 本 文 的 方法 可 以 有 效 提 高 半 高 宽 法 和 像 运 
动 法 的 视 宁 度 估算 精度 。 


4 总 结 与 展望 


本 文 分 析 了 影响 半 高 宽 法 和 像 运 动 法 估算 精度 的 各 种 因素 , 提出 了 一 种 可 以 改善 跟 踩 误 
差 ， 风 或 机 械 振动 等 因素 影响 的 方法 ， 并 使 用 实测 数据 进行 了 验证 。 我 们 有 以 下 几 点 初步 结 
论 : 

l. 使 用 半 高 宽 法 估算 视 宁 度 时 ， 使 用 更 大 口径 的 望远镜 有 利于 提高 理论 上 的 估算 精度 

2. 像 运动 法 估算 的 视 宁 度 比 谱 比 法 俩 大 ， 可 能 的 原因 是 曝光 时 间 过 长 ， 望 远 镜 口径 过 

大 ， 有 限 外 尺度 的 影响 ， 或 者 CCD 噪声 或 者 光子 噪声 的 影响 。 
3， 一 系列 斑点 图 受到 跟踪 误差 ， 风 或 者 机 械 振动 等 因素 的 影响 ， 而 表现 为 在 某 个 方向 
上 尤其 延展 开 来 时 ， 使 用 本 文 的 方法 选择 合适 的 方向 来 估算 视 宁 度 ， 可 以 有 效 改善 
O 跟踪 误差 ， 风 或 者 机 械 振动 等 因素 的 影响 。 方 向 选择 的 依据 是 : a). 此 方向 上 的 方 
差 小 ，b) . 重心 坐标 投影 在 此 方向 上 依然 较 好 的 符合 高 斯 分 布 。 

本 文 提 出 的 方法 改善 了 跟踪 误差 ， 风 或 机 械 振动 等 因素 对 半 高 宽 法 ， 像 运动 法 的 影响 ， 
提高 了 半 高 宽 法 ， 像 运动 法 和 谱 比 法 之 间 的 相关 性 。 但 是 ， 这 三 种 方法 估算 的 ro 值 之 间 仍 
然 存在 着 一 定 的 偏差 。 造 成 债 差 的 影响 因素 ， 量 化 影响 因素 与 偏差 程度 之 间 的 关系 的 方法 ， 
以 及 补偿 偏差 的 方法 是 未 来 工作 中 的 研究 对 象 。 
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Seeing estimation accuracy study based on a sequence of 


speckle images 
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University of Chinese Academy of Science, Beijing 100049 China) 


Abstract: FWHM and image motion are two frequently used methods in the field of astronomy. 
However, estimation accuracy of these two methods will be affected by tracking error, wind and 
mechanical vibration. Raising its accuracy plays a rather important role in the area of 
high-resolution imaging, site selection and seeing monitoring of observatories, etc. For FWHM 


method, we fit the long exposure image with a two-dimensional gaussian function and choose an 


appropriate direction to measure its FWHM by which we estimate the seeing parameter. A com- 
parison with spectral ratio method shows a great improvement of our approach against tracking 
error, wind and mechanical vibration. For image motion method, we employ principle component 
analysis to project center of gravity coordinates of a sequence of tracking error affected speckle 
images onto a proper direction on which we find its variance. We further utilize the variance to 
estimate the seeing parameter. A comparison with spectral ratio method also denotes a considera- 
ble improvement of our approach against tracking error, wind and mechanical vibration. 
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